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Abstract 

This report aims to present the main ideas of Regge calculus necessary to un- 
derstand the basic premise of CDT. Next, the main strategy of the CDT ap- 
proach is introduced in general terms. The main focus of this report is the 2- 
D model of CDT. The section on the 2-D model closely follows a single paper 
( [Ambj0rn fc Loll (1998)] ). While the 4-D or even 3-D case will behave very dif- 
ferently from the 2-D model, 2-D CDT can be solved exactly, and as such offers a 
better introductory exposition of CDT's methods. Higher-dimensional CDT requires 
a lot of computer simulation, and lies outside the scope of this report. 

All derivations carried out explicitly are the result of the author's independent 
work in attempting to find and prove how the results presented were obtained by 
CDT authors. Because these derivations were made explicit by the author, this 
paper can act as a guide to those who are new to CDT. 
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1 Introduction 



Physical Motivation for CDT 

One of the main goals of modern physics can be said to be the discovery of a consistent 
description of physical phenomena at all scales. Causal Dynamical Triangulation (CDT) 
is one suggested avenue of approach to this goal and thus it is of interest. 

The most attractive aspect of CDT is that it offers an approach to deriving the 
nature of spacetime from a minimal set of assumptions: the entire model arises from 
only an initial triangulation of spacetime (this puts it in the non-perturbative family of 
approaches to quantum gravity). In general, the idea of deriving what is observed from 
first principles (without having to postulate too much about the fundamental nature and 
structure of space and time) is an attractive one, and makes any model such as CDT 
that has shown initial promising results worth studying. 

Motivation for Studying CDT 

One of the main reasons for studying CDT this year is the interesting physical motivation 
presented above. Another reason for choosing CDT as a topic of research is that there is 
surprisingly little commonly known about it and there are apparently not a large amount 
of people working on it, despite it purportedly being a promising approach to studying 
how quantum mechanics and gravity may be combined. The aim was to understand why 
this may be and to understand the main ideas of the model. 

Overview and Purpose of This Report 

The goal of this MSc report and the research that shaped it is to gain a basic under- 
standing of the form and function of CDT. 

The basic idea of CDT is to discretize spacetime by approximating it with a piecewise 
linear manifold, a process called triangulation. The building blocks of the new manifold 
are simplices, which are an arbitrary-dimensional generalization of the notion of a tri- 
angle or tetrahedron. One then constructs the path integral describing the evolution of 
spacetime, which can eventually be used to obtain the physical predictions of the model. 

This report aims to present the main ideas of Regge calculus necessary to understand 
the basic premise of CDT. Next, the main strategy of the CDT approach is introduced 
in general terms. The main focus of this report is the 2-D model of CDT. The section on 
the 2-D model closely follows a single paper QAmbjOrn fc Loll (19 98)]). While the 4-D 
or even 3-D case will behave very differently from the 2-D model, 2-D CDT can be solved 
exactly, and as such offers a better introductory exposition of CDT's methods. Higher- 
dimensional CDT requires a lot of computer simulation, and lies outside the scope of 
this report. 

All derivations carried out explicitly are the result of the author's independent work 
in attempting to find and prove how the results presented were obtained by CDT authors. 
Because these derivations were made explicit by the author, this paper can act as a guide 
to those who are new to CDT. 
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2 Regge Calculus 



Before delving into CDT proper, an overview of Regge calculus is necessary in order to 
provide some of the tools required for CDT. The sources used for most of this section 
Thomson (2011)] and |Khavari (2009)] . 



were 



Regge calculus is a discretisation of general relativity in which there are no fields, just 
a triangulation of spacetime. There are two main motivations for using Regge calculus 
techniques: 

1. It offers a way of working in general relativity without using symmetries. Ordi- 
narily, significant assumptions of symmetry are required to reduce the complexity 
of relevant equations to a manageable form. 

2. It offers an approach to discretising GR, opening new avenues of approach in the 
search for a successful theory of quantum gravity. 

In the regime of Regge calculus, one considers spacetime to be made up flat (Minkowskian) 
"simplices" joined face to face, edge to edge and vertex to vertex. A simplex is a gen- 
eralization of a triangle (the 2-D simplex) to arbitrary dimension - the 3-simplex is a 
tetrahedron, a 1-simplex is a line, and a 0-simplex is a point). In this way, a smooth 
manifold can be approximated arbitrarily closely by joining any number of sufficiently 
small simplices. 

Understanding Regge calculus amounts to studying how triangulated spacetime changes 
the meaning and measure of curvature, and how that changes the form of related equa- 
tions in general relativity. The form of objects like the metric tensor under Regge calculus 
ends up being much simpler: all the information required to know the geometry is given 
by the lengths of the edges of the simplices involved. 

The simplest way to illustrate the curvature of triangulated space is to analyse the 
2-D case. Consider using a triangle with geodesic edges to probe the geometry of a 2-D 
smooth manifold. If the geometry now enclosed inside the triangle is non-Euclidean, 
the sum of the internal angles of the triangle will deviate from it. The measure of this 
deviation is the Gaussian integral curvature et, defined for a triangle t with internal 
angles a, (3 and 7: 

e t = a + (3 + j-7r. (2.0.1) 
We also need the definition of the local Gaussian curvature K at a point P: 

K(P)= lim (2.0.2) 
where At is the triangle's area. Alternatively: 

e t = [ K(P)dA. (2.0.3) 



3 




Figure 1: Triangulation of a 2-D manifold 
[Misner, Thome &: Wheeler (1973)] 

Let us now consider instead a smooth 2-D manifold whose surface we have triangu- 
lated - that is to say, we have approximated the original smooth surface by connecting 
flat triangles along their edges (see Figure 1). Using the same triangle test mentioned 
earlier to probe the geometry of what is now a polyhedral approximation to the earlier 
smooth manifold, we are led to three possible cases: 

I The triangle is placed entirely within one of the triangles composing the triangula- 
tion. 

II The triangle crosses an edge where two triangles connect, but encloses no vertices. 
Ill The triangle encloses at least one vertex. 

Since each face of the polyhedron is a flat triangle, the triangle in the first case 
encloses no curvature. The second case corresponds to taking a flat triangle and folding 
it - again no curvature is introduced. That leaves the third case: the only possibility 
of the test triangle enclosing curvature is if it encloses vertices. Thus all curvature is 
concentrated at the vertices: K(P) = if P is not a vertex (note, however, that the 
converse is not necessarily true: K{P) = can still occur if P is a vertex, as triangles 
may join at a vertex without generating any curvature). 

Of course any integral curvature is not a property of the test triangle but of the 
vertex; we will instead call it the deficiency of the vertex V, labelled ey. If a test 
triangle t encloses more than one vertex, et is the sum of the deficiencies of the enclosed 
vertices. 

An alternate way of measuring the deficiency of a vertex is by computing the sum of 
the internal angles of faces that meet at the vertex - they add up to 2tt — ey ■ This means 
that the same curvature information obtained by means of the geodesic test triangle can 
be obtained from the triangulation itself. 
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The geodesic test triangle is a useful illustrative tool in 2 dimensions, but for higher 
dimensional spaces it is much more useful to employ the notion of parallel transporting 
a vector around a closed loop. The angle cr(a) between a vector and the same vector 
after parallel transport around a closed loop a is given by: 



where A is the area of the loop o. 

In terms of our 2-D considerations above, the vector will be unchanged after parallel 
transport unless the loop contains a vertex, and upon comparing equation (2.0.4) to 
(2.0.3) we see that the angle of rotation of the parallel transported vector is exactly the 
deficit angle of the vertex (or the sum of the deficit angles of the vertices). 

The purpose of introducing parallel transport, however, was for higher-dimensional 
space considerations. In 3-D, our triangulation is formed now by 3-simplices (tetrahedra) 
connected along their flat triangular faces (2-simplices). Parallel transport of a vector 
around a loop confined to a single simplex will not change the vector; nor will the vector 
be changed for a loop that passes through a second tetrahedron but returns through the 
same face, enclosing no edges. The only way the vector can be affected is if the loop 
encloses one of the 1-simplex edges ('bones', in Regge's terminology) shared by faces of 
multiple 3-simplices. The deficiency of the edge is found by computing the sum of the 
angles between the 2-simplices meeting at the 1-simplex: the sum will be 2ir — e, where 
e is the deficiency of the edge. 

One can see how a this regime can be extrapolated to any dimension: a simplicial 
spacetime is constructed by connecting d-simplices along flat (d-l)-sub-simplices. The 
curvature is concentrated to (d-2)-sub-simplices (bones) and is computed as a deficiency 
angle by summing up the angles between the (d-l)-simplices that meet at the bone. 

Regge Riemann Curvature Tensor 

The next step is to translate some general relativity concepts into Regge calculus terms. 
The first is the Riemann curvature tensor, which will be used to obtain the Regge action. 

Each bone is a (d-2)-sub-simplex; therefore there exists a 2-D plane perpendicular 
to it. For illustrative purposes we will confine ourselves to three dimensions. Following 
[Thomson (2011)] , we consider a bundle of parallel bones (now lines) with a high density 
of bones per unit area p, each with the same deficit angle e. 

U is a unit vector parallel to the bones, and A is parallel transported around a small 
loop of area S with unit normal n and area S = nS. Decompose A into All parallel to U 
and A 1 - lying in the plane perpendicular to U . Under parallel transport around S only 
A 1 - is affected, rotated by an angle a = Ne, where N is the number of bones piercing 
the loop. In this manner, A is rotated about U by an angle a. 

From Figure 2 (6 is the angle from U to A), \5A\ = cr\A\ sinO = a\U\\a\ sin# (because 
U is a unit vector). Thus: 




(2.0.4) 




(2.0.5) 
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|A|sin0 



Figure 2: The vector A is rotated about U by an angle a. 
Thomson (2011)] 



Now a = Ne and the number of bones crossing £ is given by N = pU • £, we can 
write the above as: 



8 A = pe (U ■ £ [U x A 



or, in component form: 

5 Ay. = pee l , I/a U v A a WS 7 , (2.0.6) 

where ey Va is the anti-symmetric tensor. 

When a vector A is parallel transported around an infinitesimal parallelogram with 
sides dx and dy then <5^4 is given by: 



5 Ay = R° a pA a dx a dyP = R^A^- (dx a dy? - dx^ dy' 



(2.0.7) 



where the symmetry property of the Riemann tensor was used. The loop area £ of the 
parallelogram is given by: 



£ = dx x dy 



or, in component form: 



£u = env P dx u dy p . 



(2.0. 
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Using the identity of the anti-symmetric tensor e al3tJ, e^ U p = 5*5p — S^Sp*, 

e afS ^p = dx a dy 13 - dx?dy a . 



Substituting into (2.0.7): 



2 

Comparing this expression to (2.0.6), we find: 



s\ = he^e^K. (2.o.9) 



=> = 2peep Ua U v W. 

Multiplying both sides by e lpri and using again the identity e a ^e JPV = 5*6% - 5*8$: 

Rcrpprj Rcrfirip = 2/56 ( — U p a ) U prq , 

where U prj = e pri -yU"' '. 
Therefore we have: 

Ranpr, = peU^Upr,. (2.0.10) 

The Regge Action 

To obtain the Regge action, we begin with the regular Einstein-Hilbert action: 

Seh [g] = J d d x^~g{R - 2A). (2.0.11) 

Let us first consider only the -R-dependent part of the action. In our simplicial space- 
times the curvature is restricted to the bones; therefore, we expect that the equivalent of 
scalar curvature R under the Regge regime will vanish everywhere except on the bones 
b n . Thus we expect the independent part of the action to take the form: 

S Regge R = y~]F(b n ), 

n 

where the sum is over all bones and the function F is related to curvature and must be 
the same for all bones. We can expect F(b n ) to be proportional to the volume of the 
bone k n since the bones are homogeneous. As was stressed earlier, the only curvature- 
related variable we have is the deficit angle; we can thus conclude that F(b n ) must have 
the form: 

F(b n ) = k n f(e n ). 
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To determine the function /, consider the following argument: any bone b n can be 
represented as the superposition of two bones b n , and b n „, of the same shape and area 
such that e n = e n/ + e nff . This requires: 

/(en) = /(en/) + f{e nll ), 
which implies /(e) = Ce for some constant C. Our (independent) Regge action is then: 

SRegge R = C^2e n k n . (2.0.12) 

n 

We determine the value of C by comparing this result to the Einstein-Hilbert action 
at high bone density. Under that constraint the relationship between R and the deficit 
angle is described by (2.0.10). Using the identity U lia U tia ' = 2, we find: 

R = 2pe. 

Substituting this into the independent part of the Einstein-Hilbert action: 



L^Jd^gR^A^Jd^gpe. 



16vrG 

Recalling that p is the density of bones per unit area, corresponding to a 2-D plane 
perpendicular to a cluster of (d-2)-dimensional bones, we conclude that integrating over 
all space gives the total volume of bones: 



d d x^f^gpe = ^2k n e n , 

n 

where, as before, k n is the (d-2)-dimensional volume of the bone. Thus, 
/ d d xJ^qR ->■ — - e n k n C = — — 

J n 

by comparing with (2.0.12). Therefore: 

SRegge R = K ^ ] e ra fc ra , K = • (2.0.13) 

n 

Consider now the term proportional to A in the Einstein-Hilbert action: 

-2A f , , -2A 



16ttG 



where we have recognized that the integral gives us the volume of the entire simplicial 
spacetime, which we can express as a sum over the volume of all n-simplices composing 
the spacetime. Thus: 
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Finally, combining (2.0.13) and (2.0.14): 



SlUgge = ^2(Ke n k n -XV n ). (2.0.15) 

n 

This is the general expression for the Regge action in any simplicial spacetime di- 
mension. Some refinements which will prove useful can be made in the 2-D case, which 
we will be focusing on later on. 

Recalling equation (2.0.3), we can write: 

SRegge = £ U [ K(P)dA - XV n ] , 
n \ Jt / 

where the k n term has been dropped because a bone in a 2-D simplicial spacetime is 
a point, and therefore has no volume content. Employing the discrete Gauss-Bonnet 
theorem: 

V / K{P)dA = 2tt X , (2.0.16) 

where % is the Euler characteristic and can be written as \ = 2 — 2#, where g is the 
genus (the number of holes) of the surface. 
Thus: 

SRegge 2D = K 2lTX ~ ^ W n 
n 

=> S Re gge 2D = * 2tT X ~ N(T)XA. (2.0.17) 

In the final equation we have eliminated the sum by expressing J2 n = N(T)A, 
where N(T) is the number of triangles and A the area of each triangle. This is because 
all the triangles have the same area (in two dimensions the "volume" of the simplices 
becomes the area of the triangles) . 



9 



3 CDT Overview 



The aim of this section is to provide a general overview of the goals and methods of 
CDT, following mostly |Loll, Ambjorn Jurkiewicz (2006)| . 

In the quest for understanding of the structure of spacetime at smaller and smaller 
scales, one important conclusion we can draw is based on two premises: 

I Due to the Heisenberg uncertainty principle, probing at shorter distances introduces 
larger energy-momentum fluctuations. 

II According to general relativity, the presence of energy fluctuations will deform the 
geometry of the spacetime, imparting curvature. 

The combination of the two leads us to conclude that the structure of spacetime at 
the Planck scale must be highly curved and dynamical. 

Research into quantum gravity can be divided into two general categories: 

(a) non-perturbative: quantize the gravitational degrees of freedom, without adding any 
additional structure. 

(b) string-theoretic: the quantization of gravity appears as a by-product of a unified 
higher-dimensional, supersymmetric theory. 

CDT belongs in the first category. As will be summarized in this section, CDT 
has produced a few appealing tangible results so far: mainly, there is evidence that the 
theory has a good classical limit and the theory has provided first indications of what 
description it offers of quantum structure at Planck scale. 

The two main tools required to do CDT are path integrals and Regge calculus. The 
general form for a path integral describing the sum over the virtual paths taken by a 
particle between initial and final points x; and Xf is: 

G( Xi! x f ;t)= £ e iSP ^\ (3.0.1) 

pathsix;— >Xf 

where S paxt describes the action associated with the particle. 

For gravity, the path integral would be a superposition of all virtual "paths" the 
universe can follow as time unfolds. In this case the paths are the different configurations 
of the metric field variables g^ v {x) (accounting only for gravitational degrees of freedom; 
which is to say, ignoring matter fields). Thus we express the path integral for gravity in 
the generic form: 

G(B,gr;t)= £ e^fo-fc')], (3.0.2) 

spacetimes:gj— >gf 

where 5* grav denotes the gravitational action associated with a metric g^ v with initial and 
final boundary condition gj and gi separated by time t. The full dynamics of the system 
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can be obtained, as in the particle case, by evaluating suitable quantum operators on 
the ensemble of geometries that contribute to the path integral. CDT describes how to 
compute this path integral, and how to choose the class of virtual paths. It also provides 
technical tools for extracting physical information about the quantum geometry. 

As for the Regge calculus aspect, CDT uses the method in a manner different from 
its more classical applications. Rather than employ Regge calculus to approximate a 
spacetime, the aim of CDT is to approximate the path integral (3.0.2) as closely as 
possible, or rather to define it. 

Note that CDT does not assume that the universe is constructed from 4-simplices; 
this Regge structure is only introduced in order to make the path integral manageable. 
The edge lengths are all fixed to a common value a, and the limit a — > is eventually 
applied to the path integral to obtain its continuum limit, which is taken to be the actual 
physical prediction for the path integral of spacetime by the model. In order to achieve 
this continuum limit, the path integral is regularized: appropriate cut-off parameters 
for the configurations contributing to the path integral are introduced to make the path 
integral finite. 

The key question becomes: how do we choose which regularized triangulated geome- 
tries to select in computing the path integral? 

Previous approaches using Euclidean "spacetimes" (where, of course, all dimensions 
are treated as spatial, with no distinction being given to a "time" direction) have failed; 
they result in a Hausdorff dimension of either 2 or oo (for a 4-dimensional Euclidean 
"spacetime"), when the expected result would be 4 (however, according to some recent 
work ( [Laiho fc Coumbe (2011)] ), Euclidean DT may be able to solve this problem by 
means of an additional parameter). 

(A note on Hausdorff dimension: this is a measure obtained by comparing the typical 
linear size r of a convex subspace of a given space (e.g. its diameter) with its volume 
V(r). If the leading behaviour is V(r) ~ r du , the space is said to have Hausdorff 
dimension dn)- 

The requirement, then, is to find a path integral which allows for large short-scale 
fluctuations in curvature (accounting for the behaviour we expect to see at Planck scales, 
outlined at the start of this section), but in such a way that the resulting large-scale 
geometry does not degenerate completely (leading to a good classical limit). The ap- 
parent success of CDT in supplying this path integral lies in the imposition of causal 
rules on the building blocks; in effect, taking seriously the fact that a real spacetime is 
Lorentzian, not Euclidean, and as such contains some causal elements. 

The nature of these causal rules is as such: each spacetime appearing in the sum 
over geometries "should be a geometric object which can be obtained by evolving a purely 
spatial geometry in time, in such a way that its spatial topology is unchanged as a function 
of time" . Considerations of allowing some topology changes in time (obeying certain 
restrictions) have been made by the CDT authors, and an attempt to consider them will 
be made in the final section of this report. 
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4 The 2-D CDT Path Integral 



Now that the general strategy and goals of CDT have been presented, along with nec- 
essary tools obtained from Regge calculus, we can examine an application of the model. 
In general, obtaining results in CDT involves computer simulations and a large number 
of computations. The two-dimensional case, however, can be solved exactly, and offers 
some insight into the form taken by CDT results. The focus of the rest of this report 
is on the 2-D CDT path integral, following |Ambj0rn Sz Loll (1998)| very closely. First, 
the general structure of the 2-D triangulated spacetime is examined, and from there 
the full discrete path integral is constructed. Afterwards the coupling constants of the 
model are renormalized, and the continuum limit of the path integral is taken. Some 
interpretation of physical results is considered, after which some arguments involving 
topology changes are presented. 

4.1 Discrete Case 

The 2-D discrete CDT spacetime will have the form of closed 1-D spatial loops connected 
by triangles (Figure 3). The only geometric aspect of each spatial slice is its length, which 
is quantized in units of lattice spacing o, so that L = la, where / is an integer. Thus we 
will define spatial slices (in the discrete path integral) by I vertices (or, equivalently, I 
links connecting them). The 2-D geometry is formed by evolving an initial spatial loop 
in discrete time-steps, forming triangles. 

The rules of propagation, as expressed in the source, are as follows: each vertex i at 
time t is connected to ki vertices at time t + 1, ki > 1, by links with assigned length 
ia. To understand how the 2-D spacetime is formed, it helps to consider the process of 
evolving a spacetime in a step- by-step format: 

• Choose number of points on entrance loop: k. 

• Choose number of points on exit loop: I. 

• Each point i on the entrance is connected to ki > 1 exit points. 

• Each set of these ki points forms (ki — 1) lines on the exit loop 

=4> (ki — 1) triangles with their tip on the entrance loop per entrance point i. 

• Total number of triangles: 

— ^2i = i(ki — l) = l formed with their tip on the entrance loop (there are always 
I triangles formed with their tip on the entrance loop) 

— k triangles formed with their base on the entrance loop (formed between the 
k points on the entrance loop) 

=4> Total triangles = Yli=ii^i — 1) + k = Y^i=i h - k + k = Yli=i h = k + l 
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Figure 3: The propagation of a spatial slice from step t to step t + 1. The ends of the strip 
should be joined to form a band with topology 51 x [0, 1]. 
Ambj0rn fc Loll (1998)] 



• Finally, determine the number of ways to distribute the triangles for a given k and 
I (that is, the number of possible configurations that lead to different geometries). 
The way to calculate this number is presented a little later. 

So, the basic building blocks on this model are flat triangles with one space-like edge 
(length a) and two time-like edges (length ia). To understand where to go next, recall 
the final equation of the Regge calculus section: 

SRegge 2D =^ X ~ N(T)XA. (2.0.17) 

If we fix the topology (x = 2 — 1g is constant), the first term contributes only a 
constant phase to the path integral; we therefore ignore it and focus only on the second 
term. Computing the area of the triangles (from the general formula for the area of an 
equilateral triangle): 

. 1 2 f(iaf 1 1 2 f~ F 1 2 . [5 . 2 V5 

A = -a z \ ^ = -a\ -l - - = -a z i\ - = ia z — . 

2 V a 42V 42V4 4 

Plugging into (2.0.17): 

S = -N{T)\ia 2 ^- = Xa 2 N(T), 
where we have absorbed some constants into A. Thus: 

e iS CD T = e i\a 2 N(T)^ (4.1.1) 

The object of this section is to compute the path integral amplitude for propagating 
from a geometry of length l\ to one of length 1%. This family of expressions will have an 
iterative form, and can be fully described by three equations: 

oo 

G^(h,l 2 ) = J2G { x\h,l2;t), (4.1.2) 
t=l 



13 



oo 

G^(h,h; t) = J2 G ( x\h, I; 1) I G%\l, l 2 ;t - 1), (4.1.3) 
i=i 

Gi 1) (/ 1 ,/ 2 ;l) = ^ £ e iXa2 ^ ki - (4.1-4) 
{fci,...,fc ;i } 

The first equation is the amplitude for an arbitrary time separation of slices, derived 
by summing the second equation over t. The second equation reveals the iterative form 
the amplitude for a given t-value takes, and the third and final equation gives the base 
t = 1 case. Note that the third expression has the form anticipated in (4.1.1), recalling 
that N(T) represented the number of triangles in the spacetime, and Yli=i h = k + I 
indeed was shown earlier to represent the total number of triangles. 

It is important here to make a note about vertex marking. The above expressions 
give the physical amplitude for the cases explained in the previous paragraph. Marking 
a vertex on a loop corresponds to multiplication by the loop length factor I. This is 
because the marking of a loop means that cyclic permutations of identical triangulations 
are included in the summation (one may think of this as looking at the same configuration 
from each of I vertices on the marked loop). In (4.1.4), the curly-bracket sum includes 
cyclic permutations of the entrance loop (which has l\ vertices), and thus a factor of l\ 
appears for each unique way of triangulating the spacetime (it will be shown later, when 
the combinatoric factor J\f(k, I) is calculated, that the curly-bracket sum indeed includes 
cyclic permutations). These permutations, however, do not alter the geometry of the 
configuration, and thus do not make any physical contribution. It is for this reason 
that the l/h factor appears in the third equation: it is there to eliminate the cyclic 
permutations and give the physical amplitude. 

Throughout the rest of the 2-D CDT path integral derivation in the source (in both 
the discrete and continuous parts), vertices on loops will be occasionally marked and 
unmarked in order to simplify calculations with the understanding that the end results 
can usually be modified to give any state of marking by multiplication or division by 
Is (or Ls in the continuum limit). The form predominantly used (and which will be 
adopted for this report) is: 

G A (h,h;t) = hG^(h,l 2 -t), (4.1.5) 

where a vertex on the entrance loop has been marked to get rid of the factors of I and 
l/h in (4.1.3) and (4.1.4). 

G\(h, l 2 ; t) plays the role of a transfer matrix, obeying the properties: 

G x (h,h; h + t 2 ) = J2 h)G x (i, l 2 -t 2 ), (4.1.6) 
i 

G x (h,h; t + 1) = G ^ l > ^(J, ^ *)■ (4-1.7) 
i 

For convenience, the generating function for the Gx(l±, l 2 ; t)'s is introduced: 
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G x (x,y;t) = J2x k y l G*(k,l;t). (4.1.8) 
k,l 

We can use the above formula to rewrite (4.1.6) in terms of x and y: 

/dz 

— G A (x,z- 1 ;t 1 )G A (z,y;t 2 ), (4.1.9) 

where the contour is chosen not to include the singularities of Gx(z,y;t 2 )- 



Proof that: 



/dz 

—G x (x,z~ 1 ;t 1 )G x (z,y;t 2 ) (4.1.9) 

corresponds to the combination of the property: 

G x {k, I; h + i 2 ) = ^ G x (k, m; h)G x (m, I; t 2 ) (4.1.6) 

m 

with the generating function: 

Gx(x,y;t) = J2^ k y l Gx(k,l;t). (4.1.8) 



k,i 



Replacing the two Gxs in (4.1.9) with the form given by (4.1.8) we 
have: 



/i / -i \ rn 

fc^ZX ("J G A (fe,m;t 1 )^zVGA(n,;;t2) 

k,m n,l 



dz 

' z n - m )Gx(k,m;t 1 )Gx(n,l;t 2 ). 



=EE^v(/ 2to 

k,m n,l 

Using a change of coordinates to the effect z = re 1 ^: 

-2w ■ A, 



1 /" 2 

r n-m / 

2vr 7o 



2vr 

1 



»2tt 



27T(5 nm 
Z7T 
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=► G\(x, y; h+t 2 ) = ^2Y1 x k y l S nm G x (k, m; h)G x (n, I; t 2 ) 

k,m n,l 

= x k y l G x (k,m;t 1 )G x (m,l;t 2 ) 

k,l,m 

= E X V \^G x (k,m-M)Gx{m,l-t 2 )\ ■ 

k,l \ m J 

Recognizing the term in brackets as the R.H.S. of (4.1.6), we conclude: 
G x (x, y; h + t 2 ) = x k y l G x (k, I; h + t 2 ). 

k,l 

This is precisely the form of (4.1.8) for t = t± + t 2 . Q.E.D. 
To obtain the path integral: 

G x (x,y;l) = Y,x k y l Gx(k,l;l) 

k,l 

k,l {ki,...,k k } 

k,l {k u ...,k k } V i=l J 

=> G(x, y; g; 1) = G x (x, y; 1) = ^(^(^A/^, I). (4.1.10) 

k,l 

The summation over {k±, k 2 } indicates a summation over all possible distributions 
of values of k± , . . . , k 2 . Since each of these distributions will contribute a factor of g k+l 
(the total number of triangles is restricted to k + I at all times), we can express the 
summation as a multiplication by M(k,l), which is a combinatoric factor counting all 
possible distributions as described. 

Based on what was discussed earlier, we expect Af(k, I) to give the number of ways 
of distributing k + 1 triangles between the points of two consecutive spatial slices (Figure 
3 illustrates this well) such that no point is left unconnected (recall ki > 1). 

It turns out that this imposition (that no point is left unconnected) means that we 
can assign k triangles with their tip on the exit loop and just find out how to assign 
the remaining I triangles (which have their tip on the entrance loop). Allowing for 
maximal freedom in distributing the I triangles corresponds to there only being one 
way of distributing the k triangles (changing the placement of one of the k triangles 
always corresponds to changing the placement of at least one of the I triangles, and this 
possibility is accounted for in the various configurations of the I triangles) . 
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Figure 4: Distributing indistinguishable balls into labelled buckets. 

The problem then is finding how many ways there are to distribute I triangles with 
their tips somewhere in the set of k entrance loop points. This translates into the 
standard combinatorial problem of distributing I indistinguishable balls into k labelled 
buckets. The situation is illustrated in Figure 4. A simpler way is to think of the 
situation as a problem of how many ways there are to order a string of (I + k — 1) 
objects, where I objects are of one type and (A; — 1) objects are of another type. The 
(A; — 1) objects are the separators formed by the k buckets, and they are understood 
to be indistinguishable (it does not matter which border is first or second or third, but 
instead how many of the / balls are to the left and right of it). 

Total ways to permute the string: 

representing the fact that there are (l + k—1) items being shuffled amongst themselves, 
but / of them are identical to each other and the other (k — 1) are also identical to each 
other (so we must divide by the number of ways to order each of them to arrive at the 
true number of unique configurations) . 

Note: The ftf(k, I) formula includes cyclic permutations of the entrance loop points 
(or, in the terms of our analogy, the buckets). Since these do not produce distinct 
geometries, our final expression for G^\k,l; 1) needs a factor of 1/k to correct that, as 
is reflected in (4.1.4). Note that \M{k,l) = which is symmetric in (k, I). This 

is the overall factor that appears in G^\k, I; 1), which ( gued earlier) corresponds to 
the physical propagator for our discrete model. The propagator we work with is actually 
(4.1.5), which is not symmetric in (k,l) and is not the physical path integral (but one 
may always retrieve the physical one by dividing by k). 

Thus, returning to (4.1.10) and making the sums explicit, we have: 
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G(x,y;g;l) = 



k=0 



(i -#y) fc 



Proof that: 



E(*)' 



1=0 



{l + k-l)\ _ 1 
Z!(fc- 1)! = 



(4.1.12) 



Beginning with the R.H.S., observe that ^} z ^ k can be written as the 
(k — l)th derivative of multiplied by some factor: 



d 1 



d 1 



dzl-z (l-^) 2 ' dz{l-z) 2 {l-zf 



d 



k-l 



1 



dz k 1 \ 1 — z 
1 1 



(fc-1)! 
(l-z) k 



d: 



k-l 



1 



(l-z) fc (fc - 1)! dz*" 1 V 1 -- 2 
Now, expanding the geometric series: 



(A) 



d k 



k-l 



1 



d: 



k-l 



i=0 



dz k 1 \ 1 — z) dz k 1 
For some % = n > k — 1 (for n < k — 1 the term vanishes): 

d^ 1 



z n = Cz n-k+l 



(B) 



where C has the form: 



C = n{n- l)(n - 2)...(n - fc + 2) 

_ n(n - l)(n - 2)...(n - k + 2)(n - fc + l)(n - fc)...2 • 1 
~ (ra-fc + l)(n-A;)...2-l 



(n-ife + l)!" 
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Define I = n — k + 1 (this is our new counter with limits — > oo in the 
final, differentiated result). 



C = 



(Z + fc-1)! 



Inserting this result into (B) and (A), we have: 



1 



(1-zY 



+ k- 1)! 
U{k-1)\ ' 



Continuing with our original derivation, we have: 



G(x,y;g;l) = ^ 

oo 

-E 



fc=0 



(i -#2/) fc . 

/ oo 



v«=0 



oo / oo 



fc=0 



fc=0 



(4.1.12) 



The term with the underbrace has been added in order to exclude the contribution 
of I = 0, since that corresponds to a singularity. In this manner the also singular k = 
case also makes no contribution to our path integral. 

To get the final expression without summations: 



/ 1 \ ^ °° 

G(x,y;g;l) = ^2(gx T —) - £\ 

fc=0 v yy/ k=0 



(gx) 



l 



l - 



i - gy 



I- gx 



1 



1- gx - gy I - gx 

(i - 32/) (l - sz) - i + gx + gy 



(1 - gx)(l - gx - gy) 



G(x,y,g; l) 



(l-flfx)(l-^x-^y)" 
Combining (4.1.13) with (4.1.9), with h = 1: 



(4.1.13) 
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G x (x, y;l + {t-\)) = j> ^r- z Gx(x, z~\ l)G x (z, y-t-1) 

r dz 

= i ^z(i- 9X) (i- gx -i) G ^ y]t - 1} 

/dz cp~ 1 x 

The integrand has poles at zo = and z\ = 

To compute the integral, we use the residue theorem: 

f f(z)dz = 2m ^ Res(/; z k ) 

along with the property: 

Res^l^- F(2l) 

In our case: 



G 7 G'(zi) 



F( Z1 ) = a2xG f^ l \ G ( Zl ) = z(z(l - gx) -g) = z\\ - gx) 
l gx 

(the zo pole does not contribute because G\{zq = 0, y; t — 1) = 0). 

G' = 2z{l -gx)-g = 2g-g = g 



' F \ g 2 x \ n ( g 

Res — ;zi = G x ,y;t- 1 

G J l-gxg \l-gx 

gx n ( g n 

-Ga z ,y;t- l 



1 — gx \ 1 — gx 



G(x, y; g; t) = 9X G ( — - — , y; g;t-l 
1- gx \l- gx 



We can iterate this expression and write the solution: 



G(x,y;g;t) = F?(x)F*(x)...Fl 1 (x) 



[l-gF t ^(x)\ [1 - gF t ^(x) - gy] 1 
where F t (x) is defined iteratively by: 
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"'W' i-X.w Fo(l) = I - (4L18) 



We can verify this by examining G(x,y;g; 2): 



G(x,y;g;2) = 9% G ( — - — ,y;g; 1 
1 — gx \ 1 — gx 

gx ( g \ g 2 y 



^ V^W (l-^) (l-^-gy) 



, 2 2 

9 \ 9 xy 



1 ~ 9X) (l-T^)(l-A-^) 

7^2/ 9 2 xy 

1 {X) (l-gF 1 (x))(l-gF 1 (x)-gyY 



Continuing in this manner, 



«(.»-.,/: .9:3)= f'f(.F,. (,'))- " " ^ 



l-gxl-gx[l- 5 Fi(Fi(x))] [1 - gF^F^x)) - gy] 



= Fhx)Fi(x) l 



[l-gF 2 (x)} [1 - gF 2 (x) - gy}' 

and we can thus see that (4.1.17) holds. 

If we express (4.1.17) in terms of a fixed point, we can write it in the form: 



FAx) = F l - xF + F2t ~ 1{x - F) F= WE2 (4 119) 

t[ ) * 1 - xF + F*+i(x - FY 2g ' { 

The fixed point F is found by solving F = 

2g 

Note: The authors choose "— " sign in the fixed point formula. The justification for 
this appears to be that choosing either solution leads to no loss of generality, since both 
correspond to the same value of the physical constant g: 

9 = YTf-z- (4 - L20) 

We may verify that (4.1.17) and (4.1.19) are in correspondence by considering Ft-i(x) 
as defined by (4.1.19): 

l-xF + F 2t - 3 (x-F) 
F t -!{x) - F l _ xF + p2t _ 1 ^_ F) 
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g(l - xF + F 2t ~ 1 {x - F)) 



1 - gF t -i(x) (l-xF + F 2t ~ 1 {x - F)) - gF(l - xF + F 2t ~ 3 (x - F)) ' 

Expressing g in terms of F as per (4.1.20) and performing some straightforward 
algebra, we obtain: 
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l-xF + F a ~ 1 {x-F) 



l-gF t -i(x) " \l-xF + F 2t+1 {x- F)J ' 
We recognize the above expression as F t (x). Therefore: 

which is our iterative expression (4.1.18) (the base case Fq(x) = x can be readily verified 
to correspond to (4.1.19) by setting t = and simplifying). 

We can use this expression to express G(x,y; g;t) (as defined by (4.1.17)) in terms 
of the fixed point and obtain: 



G(x,y;g;t) 



F 2t (l - F 2 ) 2 xy 

{A t - B t x){A t -B t {x + y) + C t xy) ' 
where the time-dependent coefficients are given by: 



A t = l-F 



2t+2 



B, = F(1-F 



2/\ 



C t = F 2 (l — F 



2t-2\ 



(4.1.21) 



(4.1.22) 



Note that, by substituting in the time-dependent coefficients and performing some sim- 
plifications, (4.1.21) can be written as: 



G(x,y;g;t) 



[1-xF- F 2 *+ X (F - x)][(l - xF)(l - yF) - F 2t (F - x){F - y)\ ' 

(4.1.23) 

Expanding Fi(x)F2{x)F^{x) . . . F t -\{x) using (4.1.19): 




. . / l^{x"-F) X-^xE^F^x^F) l-xF + F 7 (x-F)" 
We can clearly see by inspection that the only surviving terms of this series will be: 

1 _ xF + F(a . _ F ) 



l-xF + FW-V+ifc-F)' 



Thus: 



F 1 (x)F 2 (x)...F t . 1 (x) = F 



t-i 



l-F 2 



I - xF + F 2t ~ l {x - F)' 



(4.1.24) 
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Now: 



1-gFt-^x) = l-gF 



l-xF + F 2t ~ 3 (x- F) 



1 - 



1-xF + F 2t ~ 1 (x- F) 
F 2 1- xF + F 2t ' 3 (x - F) ( F 
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1 + F 2 1 - xF + F 2t ~ l (x - F) 
After some algebra we find that: 

1- xF - F 2t+l {F - x) 



l-gF t _ 1 (x) = 
Similar manipulations give us: 



1 + F 2 J ' 



M 



{\ + F 2 ){l-xF + F 2t ~ 1 {x- F)) A 



(4.1.25) 



1 - gF t -i(x) - gy 



(1 - xF)(l - yF) - F (F - x)(F - y) _ N 



;i + F 2 )(l -xF + F 2t ~ 1 (x - F)) 
Substituting (4.1.24), (4.1.25) and (4.1.26) into (4.1.17): 

(1 _ i?2)2 p2 



G(x,y;g;t) = F 



2t-2 



A 



xy a2 



(1-xF + F 2t - l {x - F)) 2 (1 + F 2 ) 2 MN 



(4.1.26) 



? 2t-2 



(1-F 



2\2 



F 2 



xy 



(1 -XF; 



:iX^rF 3 f(l -xF_ 



(5^>)) 5 



F 2t (l - F 2 )'' 



xy 



MN 



Going back to (4.1.23), one recognizes M and N as being the two terms in square 
brackets in the denominator. Therefore: 



G(x,y;g;t) = 



F 2 \l - F 2 fxy 



[1-xF- F 2 *+!(F - x)][(l - xF)(l - yF) - F 2t (F - x){F - y)} 



as claimed. 

The final goal is to compute G\(h, fo; t) from G(x, y; g; t) (recall that l\, l 2 are mea- 
sures of the length of the loops; x and y have no direct physical meaning). This can be 
achieved by means of a discrete inverse Laplace transformation: 



G x (k,l;t) 



dx I dy 1 1 

n G x {x,y-t). 



Iixix J 2niy x k y 

This relation can be derived from the generating function, defined by (4.1.8): 



(4.1.27) 



G x (x,y;t) = ^2x k y l G x (k,l;t). 



k,l 
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Dividing both sides by 2irixx p and 2iriyy q : 

— — — — G x (x,y;t) = V — —x k - p y l - q G x (k,l;t). 

2mx 2my x p y q ^ 2mx 2my 

Performing a loop integral enclosing no singularities on both sides: 

-y l ~ q )G x (k,l;t). 



i — <(—--G x (x, y; t) = V ( I ^x k ~A ( I 
J 2mx J 2-KiyxPyi v y ' \J 2mx J \J 



2my* 



(4.1.28) 

Applying a coordinate change z = re 1 ^ as before to the first term in brackets: 



I dx X k- P = 1 r i+fc-p-i f 
J 2nix 2ni J 



2w 



The second term in brackets gives 5j 9 by the same method and, returning to (4.1.28) 
we have: 



j^xj^y^ Gx{X ' mt) = ^ 5kp5lq 

k,l 



G x (k,l;t) 



/dx f dy I I 
^xt^y^ ^^- 



An alternative method, used by the authors, is to rewrite (4.1.21) as a power series 
in x and y, with the final result for the discrete 2D CDT path integral being: 



rn . „ _ ~ F^B\^ m ' m( ^ ) - 1 h + l 2 -k-l (A t Q 

Ml ' 2 ' J " l 2 A\^ +2 h k\{h-k-\)\{l 2 -k-l)\\ 



(4.1.30) 
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4.2 The Continuum Limit 



Based on the similarity between our path integral formalism used in the previous section 
and the usual representation of a particle as a sum over free paths, we expect to be able 
to apply the renormalization methods of the latter in order to obtain the continuum 
limit of our discrete path integral. 

As such, we expect an additive normalization for our positive mass dimension cou- 
pling constants (i.e. the cosmological constant and the "boundary cosmological con- 
stants"): 

A = %+A, \ i = ^+X, x o = C ^ + Y, (4.2.1) 

where A, X, Y are the renormalized cosmological and boundary cosmological constants. 
Defining: 

g c = e lC \ x c = e iC ^, y c = e lC ^ , (4.2.2) 
it follows from our definitions of g, x and y in the previous section that: 

g = g c e ia2 ~ A , x = x c e iajt , y = y c e iaV . (4.2.3) 

The renormalization of the path integral, again analogously to the free particle case, 
appears as a multiplicative factor: 

G A (X, Y; T) = lim a^G(x, y; g; t). (4.2.4) 

a— >0 

The claim is that the only possible choice of r\ for which the right-hand side survives 
the limit is rj = 1. To show this, we return to a previously seen property of G\ and 
consider the form of its continuum limit: 

G x (k, I; h +t 2 ) = ^2 G x (k, m; h)G x (m, I; t 2 ). (4.2.5) 

m 

Expressing this in the shorthand form J2m f( m ) an d performing some manipulations: 
E/( m ) = I ( E/( m H = I ( E/M A/ ) dMf(M), M 

m \ m / \ m I 

Thus we have: 

G x {k, I; h+t 2 )^- a J dMG x (k, M; h)G x (M, I; t 2 ). (4.2.6) 

We require this to survive the a — > limit, therefore we need both sides to be of 
same order in a. By inspection the only way to satisfy this is if: 

G x (k,l;t) oc a, a -»• 0. (4.2.7) 



ma. 
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It will be shown later on that this is indeed the case. 

Now let us consider our expression for the generating function for G\(k,l;t): 



G x (x,y;t) = Y,x k y l Gx(k,l;t). (4.2.8) 

k,i 

Applying the process shown above but with two sums, we have: 

G x (x,y;t) ^^jj dKdLx^ ya G\(K, L; t). 



Thus we have: 

G\(x,y;t) oc \a = -, a ->■ 0. (4.2.9) 
cr a 

Therefore, in expression (4.2.4), we would need r\ = 1 so that the right-hand side is of 
order in a as a — >■ 0. 

In the paper it is argued that only for particular values of the dynamic variables will 
we obtain a non-trivial continuum limit. The authors choose \F\ 1. The reason for 
this choice seems to be getting rid of the t-dependence of G(x, y; g\ t); however, the exact 
reasons for this were not made obvious. This leads to F = e ta and thus, from our earlier 
relation between g and F: 

g c = — - — for F = e ia , aeS Usg^g c ). (4.2.10) 
y 2 cos a 1 y y ; 1 ; 

The only interesting choices for g c are stated to be g c = ±1/2, and the authors choose 
g c = 1/2 without loss of generality (explained in the next step). The reasons for this 
choice are also unclear and warrant further study. It is also determined that, given this 
choice of g c , the only interesting choices for x and y are x, y — > 1 (choosing g c = —1/2 in 
the previous step leads to the choice x, y — > — 1 and the end result is unchanged). Thus 
x c ,y c ->• 1. 

Note that the choices of values made by the authors are justified to an extent, since 
they produce a reasonable- looking continuum limit. What makes the reasons for the 
choices worth studying is the possibility of there being multiple ways of choosing lim- 
its, and the question of whether or not they result in the same continuum limit for 
G(x,y;g;t). 

In order to approach the above values from the region where G(x, y; g; t) converges, 
the renormalized coupling A is chosen to also be imaginary (A = iA). The same argument 
is followed for the so-called "boundary cosmological constants": X = iX,Y = iY . This 
leaves us with our original values g, x and y taking the form: 

9 = \e- Ka \ x = e~ Xa , y = e~ Ya . 

An important note to make is that at this stage we are talking about the Euclidean 
sector of the theory; the transition to the Lorentzian form is obtained by employing 
the analytic continuation A — > — iA, where A is our renormalized coupling constant. 
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However, the authors point out that there is a crucial distinction between this model 
and older Euclidean triangulation models: our choice of only geometries admitting a 
causal structure at the stage of our triangulation construction rules would not have 
been well justified in a purely Euclidean model (since the notion of causality implies a 
dimension of time, distinct from the rest). 

Summarizing the results of our choices of limits, we have: 

g = l -e~ Ka2 - l -ka 2 ^ , (i.e. F -> 1 - aVJ « e~ aVX ), (4.2.11) 

x = e ~ Xa » 1 - aX, y = e~ Ya » 1 - aY, (4.2.12) 

where the arrows in the first set of relations signify a redefinition of the coupling constant 
A, performed to get rid of factors of 1/2 and such in the upcoming formulae. In the 
second set we are taking the first two terms of a Taylor expansion, given that we are 
taking the a — > limit. 

The object here is substitute these expressions into (4.1.23) (the form of G(x,y; g;t) 
derived earlier) and simplify. Note that throughout the derivation the fact that we 
are taking a — > will be used to ignore next-to-leading-order terms and to express 
exponentials in a in the form given in (4.2.12). For clarity we will split G(x,y;g;t) into 
parts as such: 

G(x,y;g;t) = ^, (4.2.13) 

A = F 2 \l - F 2 ) 2 xy, B = A t - B t x, C = A t - B t (x + y) + C t xy, (4.2.14) 

where At,Bt and C% were defined earlier. Plugging (4.2.11) into the expression for A, 
we have: 

A = e -2WA( X _ e -2av / A)2( 1 _ _ aY y 

Using T = ta, we have: 

A = e- 2VXT (l - (1 - 2av / A)) 2 (l - aX)(l - aY) 

A = e- 2VXT 4a 2 A(l - aX){l - aY). 
Dropping terms next-to-leading-order in a, we have: 

A = 4Aa 2 e- 2 ^ T . (4.2.15) 

Now for B: 

B = 1 — {e- aVX ) 2t+2 - e- aVX (l - e- 2taVX )(l - aX) 
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= i _ ( e - 2 ^)( e - 2 «^) - e - aVK (l - e- 2 ^ T )(l - aX) 



Writing e -Co>/A ~ j _ c a ^/A : 



e- 2 ^ T (-1 + 2aVX + (1 - aVA)(l - aJ)) + (l - (1 - aVA)(l - aJ)) 



= e" 2 ^ T (-1 + 2a\/A + 1 - a\/X - aX + a 2 VAJr) + (l - 1 + a\/A + aJ- a 2 \/Xx) 



_ -2VAT 



(aVX -al + a 2 VXx) + (a^A + al- a 2 \/AA) . 



Dropping the a 2 terms we have: 



B = a 



(7A + x)+e- 2 ^ T (7A-x) 



Through similar manipulations, we find: 

c = a 2 [(v 7 ! + x) (v 7 ! + y) + e - 2 ^ T (VX - x) (7a - y) 

Putting A, B and C together we have: 



(4.2.16) 



(4.2.17) 



G(x,y;g;t) 
1 



4Ae - 2v /AT 



'(yx+x) + e -2VAT (yx-x^ [(yx+x) (yx+Y^j + e - 2 ^ T (yx-xj (Va-y) 

Thus, using (4.2.4) with rj = 1: 



G A (X,Y;T) = lim a G(x,y;g;t) 



^G A (X,Y;T) 



4Ae 



-2-y/AT 



(yx+x^j + e -2^AT (yx-x^ [(yx+x) (yx+Y^j + e - 2 ^ r (yx-x) (yx-yj 

(4.2.18) 

At this point several limits are contemplated (some of which will be used later on), 
and the method of obtaining Ga(Li, L2; T) forms is discussed. 
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G A (L±, L 2 ; T) can be obtained by performing an inverse Laplace transformation on 
G A (X,Y;T). Alternatively, one can take the a — > limit of G\(h,l 2 ;t) as denned in 
terms of a discrete inverse Laplace transformation in the last section. Either way we end 
up with the form: 

/ioo jy moo jy 
™ / ^ e XL, e YL 2GK{XY . T) (421Q) 
-too 2 ni J_ ioD 2m 

Consider now the T — > oo limit (this and the T — > limit will be useful for later 
derivations). Applying this limit to (4.2.18), we find: 

G A (X, Y- T) — 2 (4.2.20) 



(x + \/a) (y + \/A) 
An inverse Laplace transformation of this gives: 

G A (L U L 2 ; T) ALie ~VK(L 1+ L 2 ) e -2VKT ^ 2 21) 

In the T — > limit, applied to (4.2.18), we obtain: 

G A (X,Y;T)^ (4.2.22) 
To examine the T -> limit for G A (L ± , L 2 ;T), plug (4.2.22) into (4.2.19): 

7_ ioo 2vrz 7_ ioo 2m X + Y 

For clarity, apply the substitution X = iA', Y = iY' . 

^G A ( Ll ,L 2 ,T)^^ — £ —e e 



G A (Li,L 2 ;T) 1^ -< £ £ £ _±_ (,, 23 ) 



-oo ./— oo 

Doing the Y '-integral first: 



f°° rlV , 1 p -iX'L 2 roc iL 2 (X'+Y') 

2tt X' + Y' " 2vr ./_«, X' + Y> 



One may write the integral in the form of the exponential integral Ei{iL(X' + Y')), 
which is defined by: 

roc -t 

Ei(x) = - dt (4.2.24) 

J—x t 

which can be written in terms of the En- function with n = 1: 
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coo a — tx poo p —u 



POO g~ tx f c 

™ s L dt — = L 



du- 



ll 



so that: 



Ei(x) = -Ei(-x) 



(4.2.25) 



(4.2.26) 



To see this, consider: 



_d 
dY 



-,Ei (iL(X> + Y')) = JL { - El (-iL(X> + Y'))) = ~[ 



d roc e itL(X'+Y' 

dt 



/oo 
dt(itL) 



itL(X'+Y') 



i L J dtc ilLLK '+' > 



-iL 



p itL(X'+Y') 

iL(X' + Y') 



t=oo 



t=l 



JL(X'+Y') 

X' + Y' 



Thus: 



p -iX'L 2 r oo p iL 2 {X'+Y>) p -iX'L 2 

' dY' 6 



2vr 



-c 



X' + Y' 

JY'L 2 



27T 



00 dY' iV , T „ 1 e- iX ' L * 



2vr X' + Y< 2vr 



[£i(ii 2 (A" + rO)]~ 



Now, 



lim Ei(ix) — > in, lim Ei(ix) — > —in 

x^roo x — y — oo 



70 dY '^'L 2 1 



-iX'L 2 



2tt X' + y 2tt 



[2vri] 



/o 
-c 



rfy/ jy^2 1 _ jp -iX'L 2 

2n X' + Y> ~ 



Plugging this result into (4.2.23): 



( ; ,Ui , L, : 7 ) ^ -i I °° ^' L He-^ 



° dX' 
~2n 



jX\L^-L 2 ) 



Ga(Li,L 2 ;T)-^<S(Li-L 2 ) 



(4.2.27) 



(4.2.28) 



(4.2.29) 
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This is just as we would expect: the probability of propagating from a loop of length 
Li to a loop of length L 2 in the limit where the time for the propagation T — > should 
be zero unless L\ = L 2 , which is to say that the loop length remains unchanged. 

The authors give the final formula for G\(L\, L 2 ;T) in the general case, obtained by 
performing an inverse Laplace transform of (4.2.18): 

e^W^VK^ (4.2.30) 
V ' sinhy^T L 2 'VsinhAT )> y ' 

where h(x) is a modified Bessel function of the first kind. 

One may compute the expression for the probability of propagating from L\ to L 2 
for an arbitrary step length T by integrating the above expression over T (from — > 00). 
Thus we obtain: 



roo -VX\L 1 -L 2 \ _ -v / A(Li+L 2 ) 

G A (L 1 ,L 2 ) = / dTG A (L 1 ,L 2 ;T) = — . (4.2.31) 

Jo 2VAL 2 

One final point that must be made about the continuum limit involves the analytic 
continuation of the space-time variables. As has been mentioned, obtaining the path 
integrals for the Lorentzian theory is a matter of performing the substitution A — > —iA. 
One is naturally led to consider the analytical continuation of what we might have 
considered "time": our variable T. Attempting to use T — > —iT in (4.2.30), for example, 
yields very different (and singular) results. As the authors reason, however, this choice 
for the analytic continuation of T is wrong. 

To understand the proper way to consider analytic continuation in time, one has to 
trace the origin of the T-terms in the continuum limit equations. The term T always 
appears in the combination \J~AT , which originated by taking the continuum limit of 
F'-like terms in the discrete expressions of the previous section. 

Thus we consider instead the analytic continuation of F l . The t term is not to be 
analytically continued since it is merely a counter for the amount of iterations performed. 
It is F itself that is to be analytically continued, since its definition in terms of g (seen 
earlier) relate it to the action. To see this, first recall the definition of g: 

g = e iXa2 = e iXatai , (4.2.32) 

where we have distinguished a t , the time-direction lattice spacing, from a;, the space- 
direction spacing, for clarity. 

Looking back at (4.2.11), this means: 



9 = ~ e 
y 2 



1 - Aatai \ (\ - l -Aa t a^j , (i.e. F = 1 - = ^^^^ 



From this definition for F we can see that its analytic continuation in time would 
mean taking at — > —iaf, that is, the transformation involves converting the length of the 
time-like spacing from Euclidean to Lorentzian. This gives exactly the same result as 
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instead applying the analytic continuation to the cosmological constant, A —> —iA. Thus 
we see that thinking of T as a physical "time" -parameter is a mistake. Consequentially, as 
the authors point out, a Hamiltonian derived using it would also be physically irrelevant. 

4.3 The Differential Equation, Disk Amplitudes and more 



The aim of this section is to address a few other results discussed in Ambj0rn & Loll (1998) 



An interesting point is that one can obtain (4.2.18) by taking the continuum limit of the 
recursion relation (4.1.16). Inserting the relations (4.2.11) and (4.2.12) into (4.1.16) and 
expanding to first order in the lattice spacing one can obtain: 

^G A (X, Y;T) + -^ [(X 2 - A)G A (X, Y; T)] = 0. (4.3.1) 
The PDE is solved using (4.2.20) as a boundary condition at T = 0. The solution is: 

where X(T;X) is the solution to the characteristic equation: 

rlX 

_ = _(X 2 -A), X(T = 0)=X. (4.3.3) 

Solving this relation and plugging into (4.3.2) indeed results in (4.2.18). The authors 
then use this differential equation to obtain the Hamiltonian for the system and construct 
the solution to the "Wheeler-DeWitt equation". However, as mentioned in the article 
and at the end of the previous section, the Hamiltonian obtained in this manner cannot 
be considered physically relevant - it corresponds to the parameter T which does not 
correspond to physical time for the system. 

The next point to be made involves the so-called "disc amplitude", also called the 
Hartle-Hawking wave function, for the system. The disc amplitude gives the probability 
of a spatial slice of loop length L collapsing to a loop of zero length at arbitrary time 
(alternatively, the probability of creation of a spatial slice of loop length L from nothing 
at arbitrary time): 

W A (L) = G A (L,L 2 = 0). (4.3.4) 

We can find the exact expression for the disc amplitude in our model by applying 
the L 2 — > limit to (4.2.31). In this limit L\ > L 2 , therefore: 

e -VA\ Ll -L 2 \ _ e -^A(L 1+ L 2 ) L -y/A Ll e^ L2 - e-^ L2 _ _^ Li skWAL2 

2\/AL 2 ' 2 V / AL 2 ' VAL 2 



L ^0 ) e -VXLi 



W A (L) = e-^ L . (4.3.5) 
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The disc amplitude is then used in the source to compare to the one resulting from 
the older Euclidean model. The details of the comparison seem to draw on Dynamical 
Triangulations results too heavily for the purposes of this paper; the main point is that 
they differ. 

One more interesting result that can be drawn from our continuum path integral 
and which can be compared to the corresponding Euclidean result is the average spatial 
volume {L space ). 

In the Euclidean model, the following relation was obtained: 

G ( I u \l 1 ,L 2 ;T) oc e~^ T for T -»• 0. (4.3.6) 
We can use this to compute the average two-dimensional volume V(T): 

<F(T))< ^ W,^) ^ Gr>(Li - L2;T) 



x ~ Vat \ 4" 

=► (V(T))^ oc (4.3.7) 
For large T we expect the average spatial volume at intermediate T's to behave like: 

(Lspace)^ = ^ V -±f OC ^ (4.3.8) 

Compare now to our own model. From (4.2.21): 

G A (L 1 ,L 2 ;r)oce-^ T . (4.3.9) 
Applying the same methods as we just did for the Euclidean model, we find: 

(V(T)) oc (4.3.10) 



{L space ) oc —j=. (4.3.11) 



1 

Vt' 

This reflects the fact that the CDT quantum space-time we are working with does 
not have the anomalous fractal dimension that characterized two-dimensional Euclidean 
quantum gravity: y/A has dimension 1/ [L] , and thus our average spatial volume (which 
is of course one-dimensional in our model) has the unsurprising dimension [L\. 
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4.4 Topology 



The purpose of this final section is to mention some of the work done by the CDT authors 
in considering topology in the context of CDT. 

In |Ambj0rn k Loll (1998)| , topology change is addressed in the context of allowing 
the spatial topology to change as a function of time. This means that a "baby universe" 
is allowed to branch off from the main one at some time T, eventually disappearing into 
the vacuum - it is not allowed to rejoin the "parent" universe. This restriction (and, 
indeed, the entire consideration of this form of topology change) is imposed to permit a 
comparison with the analogous calculation in previous 2-D Euclidean calculations. The 
actual process is taken to be forbidden in CDT (at the time the paper was published) 
due to causality violations; however, some discussion of ways of possibly allowing this 
process in some form under the CDT model has been under way for some time and 
will be mentioned later on. The details of the comparison with DT and the meaning of 
many of the mathematical objects derived to allow for said comparison rely too heavily 
on knowledge of prior Euclidean work to be of relevance to this report; however, a 
qualitative consideration of the basic idea of how to represent the topology change seems 
worthwhile. 





Figure 5: A "baby universe" branches off 
locally in one time-step. 
Ambj0rn k Loll (1998)] 



Figure 6: A "baby universe" is created by a 
global pinching of the spatial loop. 
Ambj0rn k Loll (1998)| 



The process of the branching off of a "baby universe" is shown in Figure 5. Figure 
6 shows an alternative and technically simpler way to represent the process, and the 
derivation follows that representation (the continuum limit is the same in both cases). 

The modification of our approach begins with the 1-step discrete propagator, which 
takes the form: 

h-1 

G x (h, l 2 ; 1) = G$\h, h; 1) + J2 hw{ - h ~ Z ' 9)G$\h h; 1), (4.4.1) 

i=i 
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where G^ refers to the "bare" propagator without topology changes. The added term 
consists of the path integral for propagating from a spatial slice of length I to one of 
length I2 (in one step), multiplied by the discrete disc amplitude corresponding to a 
loop of length l\ — I (the length left over by the propagator). This is multiplied by h, 
representing the fact that the "pinching" may occur at any of the l\ vertices. Finally 
this term is summed over all possible ways of separating the full length l\ into two. The 
discrete disc amplitude w(l,g) (the exact nature of which is about to be discussed) gives 
the amplitude that a loop of length I collapses to a loop of one vertex (in the continuum 
limit this becomes a loop of zero length) in arbitrary time. 

The nature of w(l,g) in this formula is not straightforward; it is not the "bare" disc 
amplitude, but rather some modified version of it that accounts for topology changes 
(its exact form is unclear and beyond the scope of this section, depending heavily on 
analogues with the earlier Euclidean models). The bare disc amplitude is the discrete 
version of (4.3.5), and can be easily constructed by analogy with (4.3.4) (and (4.2.31) in 
turn for the definition of G\(Li, L2)): 



w 



(6) 



(l,g) = J2G (b) (l,h = i;g;t) = G®(!,h = !;<?)■ (4.4.2) 



Note that, as suggested above, the second argument of G^ is I2 = 1 rather than 
0, which was the argument in the continuum version of these equations; this is because 
I = is a singular case in the discrete version, which was in fact removed manually from 
the path integral in the discussion preceding (4.1.13). The contact with the continuum 
limit argument L2 = becomes evident when one considers the relation L2 = ha, with 
I2 = 1 and a — > 0. 

The purpose of the rest of |Ambj0rn Loll (19 98) | is to find an expression for Wa{X) 
under this "baby universe" regime and compare with the same result from 2-D Euclidean 
quantum gravity; this is accomplished and they are found to match after rescaling. 

Another aspect of topology and CDT involves the inclusion of a sum over topologies 
in the path integral. In [Lol l, Ambjorn Sz Jurkiewicz (2006)] , it is argued that a sum 



over topologies would cause the path integral too diverge too badly to be included in 
the expression for the path integral. |Loll &: Westra (2006)] and other papers attempt 
to find a way to circumvent this problem. This work is too involved for the scope of this 
report; the main premise is that the problem may be solvable by imposing some sort 
of causal restrictions to which geometries the sum over topologies can add to the path 
integral. 
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5 Conclusion 



In this report, the basic ideas of CDT have been introduced and its application to two- 
dimensional spacetimes presented. While two-dimensional geometry considerations are 
very different from those for three or four dimensions, the two-dimensional case is useful 
for illustrating concepts one cannot easily visualize in higher dimensional models. An 
additional reason for studying the two-dimensional case is that much of current CDT 
work (in four dimensions) relies on heavy computer simulation, but the two-dimensional 
case can be solved analytically, with its expressions and assumptions more easily iden- 
tifiable and verifiable. The details of accounting for topology changes are unfortunately 
beyond the scope of the report, but an understanding of the process in two dimensions 
is crucial to attempting to account for it in three or more dimensions. 

As mentioned in the introduction, one reason this topic was chosen for this research 
project is curiosity regarding its relative obscurity. Some of the possible reasons for 
CDT's apparent lack of popularity and low number of contributing authors may be 
its relative youth and a certain lack of marketing: there is a lack of popular and less 
mathematically intensive sources on the subject. The hope is that this report offers a 
guide to the basics of CDT and how to derive many of the expressions its 2-D application, 
and is helpful in providing a starting point from which to investigate CDT in more depth. 
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